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We introduce driven exclusion processes with internal states that serve as generic transport mod- 
els in various contexts, ranging from molecular or vehicular traffic on parallel lanes to spintronics. 
The ensuing non-equilibrium steady states are controllable by boundary as well as bulk rates. A 
striking polarization phenomenon accompanied by domain wall motion and derealization is dis- 
covered within a mesoscopic scaling. We quantify this observation within an analytic description 
providing exact phase diagrams. Our results are confirmed by stochastic simulations. 
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Understanding the physical principles governing non- 
equilibrium transport in one-dimensional (ID) systems 
has been the subject of recent interest in both bio- 
logical physics P, |3] and mesoscopic quantum systems 
Though there are fundamental differences due to 
quantum coherence effects, there is a variety of com- 
mon themes. One of them is the control of a non- 
equilibrium steady state through particle injection and 
extraction at the boundaries and coupling to some ex- 
ternal field in bulk. For example, a generic spintronic 
scheme 0] consists of sources and drains for spin injec- 
tion and extraction where the spin orientation is con- 
trolled by a tunable effective magnetic field. The analog 
in non-equilibrium statistical mechanics is the asymmet- 
ric simple exclusion process with open boundaries Q . In 
its simplest version, particles interacting only with hard- 
core repulsion move uni-directionally from the left to the 
right boundary, which are acting as sources and drains, 
respectively. These particles may either correspond to 
molecular engines like mRNA or kinesin moving actively 
along molecular tracks [B| or to macromolecules driven 
through nanoscale pores or channels [f| by some exter- 
nal field. Already this simplest conceivable model for 
collective transport exhibits phase transitions between 
different non-equilibrium steady states controlled by the 
entrance and exit rates at the boundaries It has re- 
cently been noted that a minimal model for intracellular 
transport has to account for the fact that molecular mo- 
tors may enter or leave the track not only at the bound- 
aries but also in bulk Q . Then, the bulk reservoir of 
particles acts as a gate, which can induce phase separa- 
tion [9( in a mesoscopic limit where the residence time of 
the particles is of the same order as the transit time. 

In this Letter, we introduce a generalization of the (to- 
tally) asymmetric simple exclusion process (TASEP) in 
which particles possess some discrete internal states. For 
illustration, we restrict the discussion to two states which 
are referred to as spin-up (|) and spin-down (j); see 
Fig. [TJ Inspired by Pauli's exclusion principle, we allow 
multiple occupancy of sites only if particles are in dif- 



ferent internal states. The resulting dynamics underlies 
surprisingly diverse situations. The internal states may 
correspond to distinct states of a molecular engine which 
are allowed to simultaneously occupy the same site of a 
molecular track and mutually inhibit each others motion. 
In the context of molecular or vehicular 11] traffic, 
transport on several parallel lanes may be described by 
attributing internal states to vehicles (molecular motors) 
moving on a single lane. Even further, considering hop- 
ping transport in chains of quantum dots [12] with ap- 
plied voltage, the states are directly identified with the 
spin of the electron. In this situation, the model speci- 
fied in more detail below maps to a quasi-classical version 
of a non-equilibrium Ising spin chain with nearest neigh- 
bor hopping and spin flips, where particles still respect 
Pauli's exclusion principle but phase coherence is lost. 

Motivated by this broad range of possible applications, 
we consider the following Id lattice model with L sites 
and open boundaries illustrated in Fig. [TJ At the left, 
particles with spin-up (-down) are injected at rate 
(or), respecting Pauli's exclusion principle, i.e. each site 
might at most be occupied by one spin-up and one spin- 
down state. Within bulk, particles hop unidirectionally 
to the right at rate r and may flip their spin orienta- 
tion with rate u), again under the constraint of Pauli's 
exclusion principle. Finally, the rates /3''^ control spin 
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FIG. 1: (color online) Illustration of an exclusion model with 
two internal states. Particles in states f (J.) enter with rates 
(or), move unidirectionally to the right within the lattice, 
may flip at rate w, and leave the system at rates /3 T (/3 )j 
always respecting Pauli's exclusion principle. 
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extraction at the right boundary. When interpreting our 
model as two-lane traffic [llj, the states correspond to 
lanes such that each site might at most be occupied once, 
and spin flip translates to switching between lanes. 

The quantities of interest are the state-resolved density 
p'^(x) and current profiles j''^(x) in the non-equilibrium 
steady state as functions of the spatial position x, emerg- 
ing from the interplay between external driving, coupling 
between the internal states (spin flips), and the exclusion 
principle. We identify a cooperative transition between 
a homogeneous state, where the density profiles of both 
internal states exhibit no significant spatial dependence, 
to a "polarized" state, where the density profile of one 
state changes abruptly along the track. The latter im- 
plies that the polarization p^(x) — p^(x) switches from a 
"low" to a "high" value at a well-defined position x w ; i.e., 
a domain wall forms. We find that this position can be 
tuned by changing the rates of spin injection (source) and 
extraction (drain) as well as the spin flip rate (gate). For 
the transition from a polarized to an unpolarized state 
two scenarios emerge. Either the domain wall leaves the 
system continuously through the left or right boundary 
or it exhibits a derealization transition. There are thus 
two genuinely distinct modes of switching the polarized 
state of the system on and off. 

The system's dynamics is governed by boundary (a, ft} 
and bulk (r, u>) rates, where from now on we fix the time 
scale by setting the hopping rate to unity, r = 1. Ob- 
viously, if the spin flip rate to is fast, i.e., comparable 
with the hopping rate, the spin degrees of freedom can 
be considered as relaxed such that the system's behavior 
is qualitatively the same as for TASEP without inter- 
nal states 18] . However, if the typical number of spin 



flips that a particle performs while traversing the system 
becomes comparable to the entrance and exit rates, one 
expects the interplay between these processes to yield 
interesting collective effects. In order to highlight this 
dynamic regime, we introduce the gross spin flip rate 
fl = Lu), which we choose to be of the same order as the 
boundary rates. A proper mesoscopic limit Q is defined 
by keeping Q fixed as the number of lattice sites L tends 
to infinity. 

Much of the system's behavior can already be inferred 
on the basis of symmetries and current conservation. 
Consider the motion of holes with a given internal state, 
i.e., the absence of a particle with opposite spin. These 
holes move from right to left and flip their spin state at 
the same rate u>. The dynamics exhibits a particle-hole 
symmetry: changing the notion of particles and holes 
with simultaneous interchange of "left" and "right" as 
well as a*" <-> ft'' leaves the system's dynamics invari- 
ant. In addition, there is a spin symmetry as seen by 
interchanging the spin states and the injection and ex- 



traction rates, a 



^ and 



ft . Let us now study 



the spin and particle currents passing through site i, de- 
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noted by j\ and J; = j\ + $", respectively. Since parti- 



cles are not allowed to leave or enter the system, except 
for the boundaries, the particle current Ji is strictly con- 
served and thus spatially constant, J — Ji. Unlike the 
particle current, the spin currents are not conserved indi- 
vidually. Because of spin flip processes, there is a leakage 
current from one spin state to the other. Since spin flips 
typically occur on time scales comparable to the time a 
particle needs to traverse the system (mesoscopic limit) 
this leakage current is only weak, and consequently both 
spin currents exhibit a slowly varying spatial dependence. 

Similar to TASEP, the particle current is limited ei- 
ther by the left or right boundary or the capacity of the 
bulk. The latter restricts the current to values below a 
maximal value of 1/2. For the left boundary, the current 
can not exceed Jin = cr (1 — a') + cr(l — or), while the 
right boundary constrains it to a value not larger than 
J EX = ft (I _ ft) + ft{\ _ ft), if the current is be- 
low the maximal current, it is determined by the smaller 
of the boundary currents: J = min(JiN, Jex) fli^ l- De- 
pending on which of both cases applies, we discern two 
complementary regions in the five-dimensional parame- 
ter space spanned by a*'^,ft'^ and Vl. We refer to the 
region where the total current is given by Jin as the in- 
jection dominated region (IN), while the case J = Jex 
corresponds to the extraction dominated region (EX). 
Note that both regions are connected by particle-hole 
symmetry. Since this symmetry is discrete, we expect 
discontinuous phase transitions upon crossing the border 
from the injection to the extraction dominated region, 
i.e. the IN-EX boundary. At this boundary the system 
exhibits phase coexistence, which similar to TASEP fl3l | 
manifests itself in a delocalized domain wall between a 
low (LD) and high density (HD) state of both spin states. 
Thus, based on mere symmetry arguments, we conclude 
that across the IN-EX boundary a derealization transi- 
tions appears. The simultaneous formation of a domain 
wall in the density profiles of both spin states exclusively 
occurs at the IN-EX boundary. This restriction also orig- 
inates in the conserved particle current, as can be seen 
by the following argument: For the presence of a domain 
wall in the spin-up density profile, the spin-up current in 
the vicinity of the left boundary is determined by the en- 
trance rate, jj =1 = cJ(l — cJ), while at the right bound- 
ary the exit rate specifies its value to jj =L = ft(l — ft). 
If a domain wall simultaneously forms in the spin-down 
density profile, analogous relations hold for its current. 
By conservation of the particle current, J = j\ we 
encounter the constraint 

a^(l-a r ) + a l (l-a l ) = ft(l- ft) + ft(l- ft) , (1) 

which determines the IN-EX boundary. 

Even more intriguing phase behavior emerges away 
from the phase boundary between the injection and ex- 
traction dominated regions. In particular, we find a 
broad parameter regime where a localized domain wall 
forms in the density profile of one spin state, whereas 
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the profile of the other spin state remains almost flat. 
The generic situation, as obtained from stochastic simu- 
lations and mean-field (MF) calculations (discussed be- 
low), is exemplified in Fig. [2] for the IN region. Both 
spin states enter at comparable rates, and their respec- 
tive densities approach each other due to spin flips as the 
spins traverse the system until they reach the point x w . 
There, the density of spin-up jumps to a high value (HD), 
while the one of spin-down remains at a low level (LD). 
We encounter a spontaneous polarization effect: in the 
vicinity of the right boundary, the densities of spin-up 
and -down largely differ, although they did not when en- 
tering the system. We will refer to the parameter range 
where this spin polarization effect occurs in the IN re- 
gion as the LD-HDin phase. Employing spin-symmetry 
to Fig. [5] yields a domain wall appearing in the density 
profile of the spin-down state, while one concludes from 
particle-hole symmetry that there is also a corresponding 
LD-HDex phase in the EX region. Varying the entrance 
and exit rates, one can smoothly tune the domain wall 
position as long as the IN-EX boundary is not crossed. 
If x w passes the point x w = 1 in the situation of Fig. [21 
the density of spin-up changes from the LD-HDin phase 
to the LD phase. On the other hand, x w = marks 
the transition between the LD-HDex and the HD phase. 
Regarding the domain wall position order pa- 

rameter, these transitions are continuous. 

The formation of a localized domain wall can be un- 
derstood from the continuity of the spin currents (Fig. 
[5]), which in turn arises from the only weak leakage cur- 
rent in the mesoscopic limit Q. Say the domain wall 
forms in the spin-up state. Then we have to match the 
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FIG. 2: (color online) The polarization phenomenon appear- 
ing in the IN region (q t = 0.4, = 0.2, /3 T = 0.2, (3 l = 0.45, 
and Q. — 0.5). A domain wall forms in the density profile 
of spin-up states, while spin-down stays in a LD phase. The 
spin currents, shown in the inset, are both continuous. Solid 
lines correspond to the analytical solution, while dashed lines 
indicate results from stochastic simulations on lattices with 
L = 2000 (blue or dark gray) and L = 10000 (red or light 
gray). 



currents adjacent to the left (I) and right (r) of this 
wall, which in the limit of large system size are given by 
j] /r (x) = p] /r (x)[l - p] /r (x)}. Both currents should coin- 
cide at the position of the domain wall, j/(x w ) = jl(x w ), 
while the density shows a discontinuity. Together, we ar- 
rive at the condition pl(x w ) = 1 — pj(x w ) for the domain 
wall position x w . The spatial dependence of the density 
implies that a suitable position x w indeed exists within 
a certain parameter region. 

Let us now underpin the so far general discussion by 
a quantitative analytical description. Consider therefore 
the average density p]"^ at site i for spin-up and -down 
states. The dynamical rules yield equations for their time 
evolution, which upon factorizing two-point correlations 
reduce to a closed set of difference equations; such MF ap- 
proximations have been fruitfully applied within related 
contexts 0, 



14t |15| . In a continuum limit, the differ- 



ence equations turn into differential equations, which to 
leading order in the lattice constant 1 /L take the form 

(2p T - l)d x p^ +flp l - f2p T = 0, (2a) 

{2p l - l)d x p l +f2p T - ftp 1 = 0. (2b) 

Together with appropriate boundary conditions aris- 
ing from the entrance and exit processes, they allow a 
straightforward analytic solution, the details of which 
will be presented in a forthcoming publication [l6j |. We 
have compared the analytic solution to extensive stochas- 
tic simulations through careful finite-size scaling analysis. 
Upon increasing the system size L, the densities converge 
to the analytical prediction, as exemplified in Fig. [2] The 
observed exactness of the analytical density profiles in the 
limit of large systems originates, on the one hand, in the 
exact mean-field current-density relation in the TASEP 
[4j . On the other hand, the coupling of the two internal 
states in our model locally tends to zero when the system 
size L is increased, such that correlations between them 
are washed out. The situation is somewhat analogous to 
TASEP combined with Langmuir kinetics 

The analytical approach therefore allows one to obtain 
the exact phase diagrams of the system. In particular, we 
may determine the regions where phase separation and 
thus the polarization phenomenon occurs. In general, a 
large variety of different phases appears with discontinu- 
ous as well as continuous transitions between them. As 
anticipated by our symmetry arguments, a discontinuous 
transition accompanied by a domain wall derealization 
occurs across the border between the injection and ex- 
traction dominated regions, which is given by Eq. |T]). 
Continuous transitions appear inside the IN as well as 
the EX region, when the domain wall leaves the system 
at one of its boundaries: x w — and x w — 1. The lo- 
cations in phase space where these special domain wall 
positions occur, the phase boundaries, can be obtained 
from the analytical solution of the density profiles, as the 
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FIG. 3: Phase diagram for equal injection rates cJ = or = a. 
Phases for spin-up and -down are shown in panels (a) and 
(b), respectively, as a function of a and ft for fixed values 
ft = 0.3 and f2 = 0.15. Phase separation (shaded areas) 
arises in the IN as well as the EX region. The delocalization- 
transition line (thick line) is identical in (a) and (b). Thin 
lines correspond to continuous transitions. 



latter reveals the position x w ; details will be presented 
elsewhere lid, . 

For illustration, we focus on the special case of equal 
entrance rates, or = or = a, which already exhibits 
the main features of the possible non-equilibrium steady 
states. In particular, when a domain wall occurs in the 
density profile of one of the spin states within the IN 
region, it lucidly shows the polarization phenomenon: 
both spins have equal densities in the vicinity of the left 
boundary, but strongly differ at the right one. The phase 
diagram for the spin- up and spin-down states, resulting 
from the analytical solution, is presented in Fig. [3] In 
the two-lane interpretation of the model, the states refer 
to the upper and lower lane. For both internal states the 
same first order line marks the transition between the 
IN and EX region. The polarized state (shaded area), 
where a domain wall appears either in the spin-up or - 
down state, intervenes between the LD and HD phases. 
There are continuous transitions from pure (LD, HD) to 
coexistence phases as the domain wall enters or exits the 
system at the boundary. The lines marking these tran- 
sitions intersect the IN-EX boundary in a multicritical 
point ,4, where all boundary rates equalize, a — ft — ft. 

Consider now a horizontal path through the phase di- 
agram below the multicritical point A. For small val- 
ues of the injection rate a the system is in the injection 
dominated region, and both spin states (lanes) are in a 
homogeneous LD state. Then, upon crossing the phase 
boundary for the spin-up state, the system switches to a 
polarized state, similar to Fig. [21 such that a domain wall 
appears in the density profile of the spin-up state (upper 
lane), entering continuously from the right boundary; for 
the same parameter range the spin down state stays in a 
homogeneous LD state. Approaching the IN-EX bound- 
ary this domain wall delocalizes, and upon crossing re- 
localizes again, but now as a domain wall in the spin- 



down state (lower lane) . As the spin- up state has turned 
to a HD phase, we encounter polarization near the left 
boundary Crossing the IN-EX boundary, the system 
thus switches its polarization from the right boundary 
to the left one. A further increase of the injection rate 
finally shifts the position of the domain wall to the left 
boundary such that the system remains in a HD phase 
for both spin states (lanes). For a path through the phase 
diagram above A, similar arguments hold. 

Two of the lines that mark continuous transitions 
are readily obtained. The transition from the LD to 
the LD-HDin phase is determined by x w = 1; since 
p*(x) — p^(x) — a for x < x w , it is located along the di- 
agonal, a = ft , for spin- up state and parallel to the ver- 
tical axis, a = ft, for the spin-down state. We emphasize 
that these phase boundaries in the injection dominated 
regime do not depend on the magnitude of the gross spin 
flip rate f2; i.e., qualitatively tuning the system's state is 
possible only upon changing the injection or extraction 
rates. The lines corresponding to continuous transition 
in the EX region are more involved [l6| ■ Its most notable 
feature is that the width of the polarized phase decreases 
with increasing spin flip rate fl, until it finally vanishes 
in the limit — > oo. 

The authors are grateful for insightful discussions with 
Felix von Oppen and Ulrich Schollwock. 
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